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Exchange in a silicon-based quantum dot quantum computer architecture 
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In bulk silicon, intervalley electronic interference has been shown to lead to strong oscillations 
in the exchange coupling between impurity electronic wavefunctions, posing a serious manufactura- 
bility problem for proposed quantum computers. Here we show that this problem does not arise 
in proposed architectures using Si/SiGe quantum dots because of the large in-plane strain in Si 
quantum wells together with the strong confinement potential typical of heterostructures. 
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Introduction. The proposal to build a quantum com- 
puter with qubits that are single electrons confined in 
gated silicon quantum dots in silicon/silicon-gcrmanium 
heterostructures P, 01 is attractive because of the self- 
aligning properties of quantum dots as well as the long 
coherence times of electron spins in silicon 0, |^ . How- 
ever, one complication of silicon compared to a direct 
band gap material such as gallium arsenide is the valley 
degeneracy in the band structure. As Rcfs. 0, S S 
point out, interference between valleys causes charge os- 
cillations that can cause the coupling between qubits to 
vary rapidly in space. Since the fundamental two-qubit 
gate is achieved by varying the overlap of electron wave- 
functions centered at different locations, obtaining reli- 
able gates would appear to require controlling the indi- 
vidual wavefunctions on the scale of the oscillations and 
not the much longer length scale describing the variation 
of the envelope of the wavefunctions. This sensitivity 
has severe implications for the feasibility of constructing 
a reliable quantum computer using silicon. 

Here we present theoretical arguments and atomistic 
simulations that show that the large in-plane strain 
present in a silicon quantum well, together with the 
strong potential confining electrons in the quantum well, 
combine to eliminate the problem of the fast oscillations 
in the overlap of wavefunctions of electrons centered at 
different locations in a Si quantum well in a Si/SiGe 
heterostructurc. The result depends on the strong hct- 
erostructure confinement present in the gated quantum 
dot architecture, and does not apply to other Si-based 
quantum computing schemes [ittlllLll^ . 

Effect of Strain on Electron Wavefunctions in Si 
Quantum Wells. Unlike gallium arsenide, unstrained 
silicon has six-fold degenerate conduction band minima 
located along the [100], [010], and [001] directions, about 
85% of the way to the boundary of the Brillouin zone. In- 
terference between these valleys causes fast oscillations in 



the electron density, on the scale of the interatomic spac- 
ing. In bulk silicon the oscillations occur along the a;, y, 
and z directions. These fast oscillations make controlling 
the exchange coupling accurately, which is necessary to 
make a two-qubit gate 0] , much more difficult d, Q, . 

In Si/SiGc heterostructures the large in-plane strain 
present in the quantum well reduces the six-fold valley 
degeneracy to a four-fold and a two-fold one. For a typ- 
ical quantum well in the x-y plane, perpendicular to the 
[001] crystal axes, the two lowest energy valleys arc at 
kx ^ 0, ky ^ 0, and ifc^, with kz ^ 0. For typical 
heterostructures, the minima with nonzero kx or ky have 
energies more than 0.1 eV higher than the minima with 
nonzero kz, so that at the low temperatures at which a 
quantum computer would operate, only the valleys with 
kx = ky =0 are relevant. Therefore, in the plane 

of the quantum well the charge density does not oscillate. 
However, the two degenerate valleys along ±kz do lead 
to oscillations along the z direction. Therefore, strain 
alone docs not completely remove the extreme sensitiv- 
ity of the exchange coupling to small changes in position, 
especially for impurity based qubits 0, 0- 

Below we show that in a Si/SiGe quantum dot quan- 
tum computer wavefunction oscillations along z do not 
lead to problems with controllability of the exchange in- 
teraction because the oscillations of the wavefunctions 
of different qubits arc aligned by the strong confinement 
potential of the quantum well. The argument is quite 
general, and the result is robust even in the presence 
of gate potential fluctuations and imperfections such as 
quantum well width variations. 

Role of Confinement Potential. The exchange in- 
teraction does not exhibit fast oscillations in quantum 
wells because the oscillations along z are aligned by the 
strong quantum well potential. The physical reason for 
the alignment of the charge density oscillations perpen- 
dicular to the quantum well plane is that the quan- 
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turn well confinement potential varies on a much shorter 
length scale than any in-plane potential variations. (Re- 
call that typical quantum wells have 10 nm widths and 
potential steps occurring on the single unit cell scale 
and heights ^ 0.1 eV [13, while the potentials defining 
the quantum dot vary on length scales of many tens of 
nanometers.) Because of this separation of length scales, 
it is appropriate to use the Born-Oppenheimer (B-0) ap- 
proximation |l6l IT^ . We wish to find the lowest energy 
eigenstate ipE of the time-independent Schrodinger equa- 
tion of an electron in a single quantum well, 

^2 ^ Q2 Q2 Ql ^ 1 



2ra \9a;^ dy 



= E^Pe{x,v, z) .(1) 

We write the potential as the sum of three 

contributions, the atomic potential T4(a;,j/,z), the con- 
finement potential Vc{x,y^z)^ and the residual potential 
Vr{x,y^ z)^ which accounts for potentials from external 
gates and possible weak impurity potentials (e.g., from 
dopant inhomogcneities in the modulation doping layer). 
For simplicity, initially we will ignore the atomic po- 
tential Va{x,y,z) and consider a system with confine- 
ment and residual potentials whose sum varies much 
quickly along the z direction than in the x — y plane. 
We follow the usual B-0 argument 0, 0| and write 
"fpix, y, z) = ^o{z; X, y)x{x, y), where ^o{z; x, y) describes 
the ground state wavefunction as a function of z at given 
x and y, and x(^i?/) describes the in-plane variations. 
The B-0 prescription is to first find $0 by solving 



- — — + V{x,y,z) 



and then determine x(x, y) from 



<^o{z;x,y) = Veffix,y)<^o{z;x,y), 

(2) 



x{x,y) = Ex{x,y). 



(3) 

Eq. for the z-dcpcndcncc of the wavefunction de- 
pends only on the local value of x and y. The wave- 
functions obtained for two different residual potentials 
in the x — y plane, Vriir) and Vr2{r) (corresponding to 
electrons centered on two different quantum dots), are 
^o{z;x,y)xi{x,y) and ^o{z;x,y)x2{x,y). Here, xi and 
X2 are distinct, but $0 is the same for the two wave- 
functions. Therefore, while the wavefunction can vary 
quickly as a function of z, this z-dependence is identi- 
cal for wavefunctions centered at different locations in 
the X ~ y plane. The overlap matrix element between 
electrons centered at fi and r2, (n.|^2) is 



(ri|r2)= / dxdy x*i{x,y)x2{x,y) 



dz\^o{z]x,y)\^ 



The term in brackets does not depend on the residual 
potential Vr{x^y), so it does not change when the gate 
potentials are varied. In other words, the strong quantum 
well confinement "locks" the variations perpendicular to 
the quantum well plane. Since, as discussed above, the 
wavefunctions in these strained quantum wells have no 
oscillations in the x ~ y plane, the scale of all variations 
in the exchange coupling is the quantum dot size, typ- 
ically of order at least 100 nm. Therefore, the valley 
degeneracy does not lead to additional problems in con- 
trolling the exchange interaction in a Si/SiGe quantum 
dot computer. 

To include the atomic potential Va(x, y, z), we use the 
envelope approximation Il9j and write the wavefunction 
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where Uj{r)e'^^^''^^ is the Bloch wave at the minimum of 
the j*'' valley, the aj describe the amplitude of the con- 
tribution from each valley, and the Fj (r) are "envelope" 
functions. Because band offsets in typical heterostruc- 
turcs are ^ 0.2 cV, while atomic potentials are several 
eV, the variations induced by the confinement potential 
are much slower than those induced by the atomic poten- 
tial, and therefore the envelope function itself satisfies the 
Schrodinger equation, Eq. (QJ Thus, within the en- 
velope approximation, the argument given above with no 
atomic potential applies with no modifications. 

In conclusion, we have shown that the strong confine- 
ment potential in Si/SiGe quantum wells causes the ex- 
change coupling to depend smoothly on quantum dot sep- 
aration. 

Tight-Binding Model Calculation. In explicit 
support of the aforementioned conclusions we have nu- 
merically computed the behavior of the exchange cou- 
pling as a function of the separation between quan- 
tum dots in Si/SiGc hetcrostructurcs by using the quan- 
titative nano-electronic modeling tool NEM0-3D [20j . 
NEM0-3D describes the electron Hamiltonian in the 
framework of an sp^d^s* nearest-neighbor empirical 
tight-binding model, which allows us to incorporate ex- 
plicitly the effect of well- width variation, interface rough- 
ness, and strain at the atomic level. The atomistic de- 
scription is essential to verify the aforementioned argu- 
ment that the electron wave function varies smoothly in 
the X — y plane without fast oscillations at the atomic 
scale. We use the empirical tight-binding parameters of 
Ref. [iJl , which reproduce both the band edges and effec- 
tive masses of the lowest conduction band and the highest 
valence band to within less than 5%. The effects of strain 
on the band edges and effective masses are incorporated 
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FIG. 1: Schematic cross section of the modeled system which 
consists of a Si quantum well with a width of 4.4 nm and a 
Sio.rGeo.a barrier material with a width of 2.8 nm below and 
above the Si well. The width of one half of the Si well is 
increased by one monolayer to imitate a well-width variation 
arising from growth inhomogeneity. 

into the model by modifying the tight-binding parame- 
ters with the Lowdin orthogonahzation procedure, the 
Slater-Koster table, and t he g eneralized version of the 
Harrison's scaling law |22| . 

Figure^shows the modeled system, a strained Si quan- 
tum well with a width of 4.4 nm and a relaxed Sio.yGeo.a 
barrier buffer with a width of 2.8 nm below and above 
the Si well. The supercell dimension in the x-y plane is 
40 nm. To simulate the effect of the well- width variation 
due for example to a miscut substrate, the Si well width 
is increased by one monolayer in half the supercell in the 
X direction. The external gate potential Vg{x,y,z) for 
the lateral confinement is approximated by a harmonic 
potential in the x — y plane and a linear potential along 
the z direction: Vg{x,y,z) = A{x'^ + y'^) + Bz, with A and 
B chosen to be 0.1 meV/nm^ and 20 meV/nm, respec- 
tively, in agreement with electrostatic potential calcula- 
tions AH other potentials Vc{x,y,z) and Va{x,y,z) 
are included in the tight-binding Hamiltonian via the 
couplings between the tight-binding basis orbitals. By 
diagonalizing the resulting Hamiltonian, the ground state 
electron wave function Tp{r) for a single quantum dot is 
obtained. The two-electron wave functions in a system of 
two quantum dots are prepared by superposing the single 
electron wave functions centered at each of the two quan- 
tum dots. Two types of two-electron wave functions can 
be constructed: symmetric and antisymmetric states. 

The exchange coupling J of the two-electron system 
or, equivalently, the energy difference between the sym- 
metric and antisymmetric states is given approximately 
by i 

J{R) = j dxi dyi dzi dx2 dy2 dz2 

rinWiri ~ R) .^^ ^ ijin - R) ^(ri), (6) 
qri - r2 

where R is the relative distance vector between the cen- 




Interdot Separation along [100] (nm) 



FIG. 2: Exchange coupling as a function of the interdot sep- 
aration along the [100] direction for two quantum dots that 
are electrostatically defined in a Si quantum well surrounded 
by a Sio.rGeo.s buffer. The inset shows the overlap between 
single-electron wave functions ip{r) centered at each of the two 
quantum dots, i.e. \ J dx dy dz ip* {f— X)ip{r)f , as a function 
of the interdot separation X along the [100] direction. Both 
the exchange coupling and the overlap (inset) smoothly in- 
crease as the interdot separation decreases. This result clearly 
demonstrates the absence of atomic-level oscillations. 

ters of the two quantum dots and e is the dielectric con- 
stant. This exchange integral is further expanded into 
integrals involving tight-binding basis orbitals. The de- 
tails of the expansion and the calculation of the integrals 
can be found elsewhere [23j . 

Figure |21 shows the exchange coupling as a function 
of the interdot separation along the [100] direction. As 
the interdot separation decreases, the exchange coupling 
smoothly increases without atomic-level oscillations. To 
show the origin of this behavior, we also plot the over- 
lap between single-electron wave functions centered at 
each of the two quantum dots in the inset of Fig. [21 As 
expected, the overlap smoothly increases with the de- 
crease of the interdot separation. These results clearly 
demonstrate the absence of atomic-level oscillations in 
the single-electron wave function, and consequently the 
absence of fast oscillations in the exchange coupling in 
the two-electron wave function. 

Summary. We have shown that the exchange inter- 
action between two qubits composed of quantum dots in 
a Si/SiGe heterostructure exhibits a smooth dependence 
on qubit separation. The origin of the difference between 
the behavior in quantum dots and in bulk silicon impu- 
rities is the strain and the strong confining potential in 
the heterostructure. 
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